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RESPONSE OF CONVECTIVELY CONTROLLED BURNING TO NONLINEAR DISTURBANCES 


"by Marcus F. Heidmann 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

A process where energy release rate is proportional to the square root 
of the Reynolds number of a burning particle or surface is used to demonstrate 
some distinctive dynamic properties of a disturbed convectively controlled 
burning process. The instantaneous energy release rates caused by assumed, 
periodic flow field disturbances are numerically evaluated and examined. 
Correlation coefficients which express the energy released in r. phase with the 
pressure disturbance are evaluated and analytical solutions for these 
coefficients are derived. 

The response of the process is shown to be highly sensitive to the harmonic 
distortion of the disturbance. For some disturbances the energy released in- 
phase with the pressure disturbance is an order of magnitude larger than that 
for linear (sinusoidal) disturbances. The results show that harmonic distor- 
tion increases the coupling between the burning process and the flow field. 

This amplifying effect of harmonic distortion is suppressed when a disturbance 
becomes steep-fronted; therefore, steepening can act to limit the equilibrium 
amplitudes in disturbed systems. Amplification is also suppressed when the. 
velocity disturbance is pha^e shifted with respect to the pressure disturbance 
or when the process is exposed to a high steady velocity. 
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INTRODUCTION 

The energy release process in many reactive systems depends on the con- 
vective forces on a "burning particle or surface. For such systems, "both the 
density and velocity associated with nonsteady flow will affect the rate of 
energy release. The purpose of this paper is to examine and evaluate the 
dynamic coupling "between a convectively controlled "burning process and a 
nonsteady flow field where pressure, density and gas velocity vary periodically. 
Emphasis is placed on distorted (nonsinusoidal) flow field disturbances which 
include these with "shock-like" wave shapes. 

For most reactive systems, linear (small sinusoidal perturbation) analyses 
provide much of the basic insight into the behavior of the disturbed system. 
Linear analyses are valuable because one can readily comprehend the physical 
behavior described by the analysis. However, the distortion and gas motion 
which accompany many disturbances are neglected in linear analyses. The 
practitioner in real systems should know and understand the significance of 
these deviations from linearity. The method of analysis employed in this study 
is directed toward the objective of expanding upon the insight which linear 
analyses provide by including the affects of distortion and gas motions in the 
evaluation of dynamic properties. 

A process where energy release rate is simply proportional to the 
Reynolds number of a burning particle or surface to an exponential power is 
used to demonstrate some distinctive dynamic properties of the convectively 
controlled burning process. The instantaneous energy release caused by 
assumed flow field disturbances are numerically evaluated and examined. A 
correlation coefficient which expressed the energy released in-phase with the 
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pressure disturbance is evaluated and analytical solutions for the coeffi- 
cient are derived. The effect of distortion and gas velocity on the coeffi- 
cient is discussed with regard to the growth and decay of the disturbance. 

The method of analysis and some of the results have been reported pre- 
viously (1,2) in studies of combustion instability in liquid propellant 
rocket engines. The broader implications of these previous studies are re- 
viewed without the application to a specific reaction system. Additional 
analytical solutions and evaluations are presented for steep-fronted and 
"shock-like" disturbances which are frequently encountered in a variety of 
reactive systems. 

MODEL 

The model used to evaluate the response properties of a convectively 
controlled burning process to flow field disturbances includes the following 
elements : 

Energy Release Process 

A process is assumed where energy release rate is given by: 

V - VRe)l/S . (1) 

An energy release rate dependent on Reynolds number, as expressed by 
equation (l), characterizes the convective heat and mass transfer mechanisms 
which dominate many energy release processes. The gas viscosity, p, and the 
characteristic dimension, D, are assumed constant. The velocity parameter U 
is the magnitude of the relative velocity seen by the particle or surface. 

For this study a two-dimensional flow field which is steady in one dimension 
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The angle cp n (harmonic phase angle) gives the phase relation between 
the harmonic components. As -in linear analyses, pressure and density dis- 
turbances are assumed to be in-phase with each other. The angle 0 n (velocity- 
pressure phase angle) expresses, the phase relation between the velocity and 
pressure disturbances. For traveling waves, which will often be discussed. 


u n cos (neat 
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the velocity and pressure are assumed in-phase as in linear analyses. 

Equation (4) is simplified for some specific evaluations in order to 
reduce the number of variables. The assumptions that p n = irp n and u n = cp n 
are made. These are approximations adopted from linear theory and adequately 
characterize the gas dynamics for the demonstration purposes of . this study. 

Specific evaluations are made for two types of harmonically distorted 
disturbances. In one type the distortion is limited to two harmonic compo- 
nents. In the other type multiharmonic distortion is considered where the 
harmonic coefficients are assumed ordered by p n = p^ and the phase angles 
(Ph and 0 n are assumed constant. Such highly ordered harmonic series can 
describe the "spiked" pressure wave shapes observed in strong traveling 
acoustic modes (4) as well as "shock-like" shapes. 

Response Properties 

The response properties are expressed by correlating parameters which 
relate the perturbation in energy release rate to the imposed disturbances 
in the flow field. One of the more important correlating parameters is the 
component of the energy release rate perturbation which is in-phase with the 
pressure disturbance. According to Rayleigh's criterion for heat driven 
waves, this in-phase response gages the ability of the process to amplify or 
attenuate the flow field disturbance. 

For this analysis the component of the energy release rate perturbation 
in-phase with the pressure disturbance will be called the in-phase response 
factor, ft. It can be extracted from the energy release rate perturbation 
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and normalized by a correlating procedure (l) defined by: 

n 2it 

W'P' dent 

Jo 

R = n2rt ( 5 ) 

Jo ( p ’> 2 

By linear analysis, it can be shown that energy is added to the flow 
field disturbance when ft is positive and energy lost when negative. Growth 
or decay rate of the disturbance depends on the magnitude of ft. In unstable 
rocket combustions, where energy is continually lost due to exhaust flow, 
ft must usually be greater than about unity for sustained oscillations. 
Whether such linear criterion apply directly to nonlinear disturbances has 
not been rigorously established. However, applications such as that made 
to Priem and Guentert's nonlinear studies in reference 1 show that criterion 
established by linear analyses acceptably predict nonlinear behavior. 

Analytical Solutions 

Approximate analytical solutions for the response can be derived by ex- 
panding the burning rate expression (eq. 3) in a Taylor series. Details of 
this method are given in reference 1. By neglecting harmonic content higher 
than second order in the Taylor series and cross products of harmonic coeffi- 
cients higher than the fourth order the following general solution for ft is 
obtained for cp = 0: 
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In this solution, the amount of second harmonic distortion is expressed 
as a ratio of second to first harmonic amplitude^, p /p T . The first term in 

equation (8), l/2y> is the linear response - the response to small sinusoidal 

! 

disturbances. 

The solution for R to second order in p^ for the multiharmonic dis- 
turbances where p n = p^ and cp and 0 are constants is given by: 
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In all of the solutions, the response factor £ is a maximum when the 
phase angles cp and Q are equal to zero. For equation (9) this condition 
of zero angles characterizes "spiked", traveling type disturbances. 

DISCUSSION OF RESPONSE PROPERTIES 

In the first part of this discussion the effect of the distortion in a 
disturbance on the instantaneous ' energy release rate will be described to 
illustrate some distinctive properties of the convectively controlled burning 
process. The instantaneous energy release rates were evaluated numerically 
with the aid of a digital computer. The discussion of the response factor £, 
which follows, is also supplemented with some numerical evaluations which help 
establish the precision of the approximate analytical solutions for £. 

Instantaneous Energy Release Rates 

The effect of harmonic distortion on the instantaneous energy release 
rates is best illustrated by considering the multiharmonic disturbances 
where p n = p^. For these disturbances, harmonic distortion increases with 
pressure amplitude. Figures 1-3 shows the change in the pressure wave 
shape with amplitude and compares the numerical solution for the energy re- 
lease rate with the pressure disturbance for several conditions. 

Figure 1 illustrates a type of reference condition in which velocity 
disturbances are neglected. The energy release rates are equivalent to 
those for a process which is only density sensitive (w ~ p-*-/^). Without 
velocity oscillations, the energy release rate essentially follows the pres- 
sure disturbance but at a reduced amplitude. The peak energy release rate 
varies from about one -half the peak pressure for disturbances which are 
sinusoidal (fig. 1(a)) to much less than one-half for the "spiked" distur- 
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bances (fig. 1(d)). There appears to he a loss in response with an increase 
in the harmonic distortion for a process which is only density sensitive. 

Figure 2 illustrates the importance of velocity disturbance for the 
same pressure disturbances considered in figure 1. The effects of velocity 
disturbances are very pronounced for sinusoidal type disturbances (figs. 2(a) 
and (b)). The energy release rate perturbation is much larger than the 
pressure disturbance but exhibits two peaks. Two peaks exits because the 
convectively controlled process is sensitive to the magnitude and not the 
direction of the velocity vector and the magnitude is high at both high and 
low pressure conditions. With regard to coupling between the flow field 
and the process, the .two peaks appear to have counteracting affects, .that 
is, energy is released first in-phase and then out-of -phase with the pres-, 
sure. 

The correspondence between energy release rate and the pressure dis- 
turbances in figure 2 progressively improves with an increase in pressure 
amplitudes, When the amplitude and distortion are high (fig. 2(d)) the 
'Ispiked" disturbances cause a very similar disturbance in release rate. 

The second peak in release rate seen at low amplitudes is suppressed but. 
the peak which is in-phase with the pressure peak remains large. The energy 
release rate perturbation appears as an amplified pressure disturbance. 

Figure 3 illustrates the dynamic behavior for disturbances which be- 
come steep-fronted and "shock-like" with an increase in pressure amplitude. 
The results are very similar to those for the symmetrical type "spiked" 
disturbances illustrated in figure 2. An increase iri amplitude and harmonic 
distortion improves the correspondence between the energy release rate and. 


pressure disturbance. 
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These effects of harmonic distortion on the dynamic behavior of con- 
vectively controlled burning can be discussed from a more quantitative stand- 
point by a comparison of the in-phase response factors which relate the en- 
ergy release rates to the pressure disturbances. 

Response Factor 

The effect of an increase in the harmonic content of a disturbance on 
the response factor ft is the best illustrated by initially considering 
disturbances which contain only first and second harmonic components. The 
analytical solution given by equation (8) applies to such disturbances. 

Figure 4 shows the effect of second harmonic distortion on the wave shape 
and on. the in-phase response factor ft for disturbances where the harmonic 
phase angles <p and 6 are equal to zero. The response factor exhibits 
optimum properties as a function of both pressure amplitude the root- 

mean-square amplitude of the disturbance) and second harmonic content 
(Pg/Pp* the ratio of second to first harmonic pressure amplitudes). Maximum re- 
sponse occurs' with' a second harmonic Content of about 0.8 at an amplitude of 0.02. 
This maximum is about an order of magnitude larger than that for linear 
(sinusoidal) disturbances which is shown by the p /p = 0 curve in figure 4. 

The result implies that second harmonic distortion can have a very large 
effect on the growth of a disturbance. 

The effect of increasing the harmonic content in disturbances with 
multiharmonic distortion is shown in figure 5. The analytical solution 
given by equation (9) applied to such disturbances. For figure 5, the har- 
monic and velocity-pressure phase angles are equal to zero. These distur- 
bances are of the type considered in figure 2. Figure 5 shows that an in- 
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crease in pressure amplitude which increases the harmonic distortion causes 
an increase in the response factor for any constant value of relative 
Mach number, (Av/c). At low amplitudes, the response is the linear value 
of l/2y and it asymptotically approaches a value of about 3 times this 
linear value at high amplitudes. The result implies that in some systems 
a disturbance above some finite amplitude may be needed before the response 
is sufficiently large to cause a disturbance to grow ih amplitude. Rocket 
engine combustors exhibit such .behavior. 

Figure 5 also shows the effect of the steady velocity (Av/c) on the 
response. The convectively controlled process is more sensitive to velocity 
oscillations when the steady velocity environment of a particle or surface 
is small. Figure 5 shows that the response factor increases with a decrease 
in the steady velocity. For any real process, however, there is a lower 
limit . to this steady velocity where convective forces no longer control the 
rate of energy release and this analysis is no longer applicable. 

Analytical and numerical solutions for the response factor can be com- 
pared in figures 4 and 5. Such comparisons have shown that the analytical 
solutions generally over-predict the response by about 10 percent but ad- 
equately express the functional dependence of the response factor on the 
distortion variables. 

The harmonic phase angle cp affects the response factor as shown in 
figure 6, As cp increases, the disturbance changes from being symmetrical 
about the peak as shown in figure 2 and becomes steep-fronted like those 
shown in figure 3. Such wave steepening caused by an increase in cp is 
shown in figure 5 to suppress the increase in £ with an increase in pres- 



12 


sure amplitude. At an angle of jt/2 the response remains near the linear 
value of l/2y even at high pressure amplitudes. This result implies that 
wave steepening may he a mechanism whereby pressure disturbances reach 
limiting amplitudes in reactive systems. That is, if a disturbance steepens 
and becomes "shock-like" as it grows in amplitude the energy added in-phase 
with the disturbance would decrease and eventually reach a condition where 
no more growth is possible. 

For all the examples discussed thus far, the velocity disturbances 
have been in-phase with the pressure disturbance (e = 0) . Figure 7 shows 
the effect on the response factor of varying this phase relationship. An 
increase in the velocity-pressure phase angle 6 is shown to suppress the 
in-phase response factor in a manner similar to that for the harmonic phase 
angle cp, but the effect is less severe. At an angle of n/2, the response 
ft remains above the linear value' of l/2y at high amplitudes. In reactive 
systems, the velocity-pressure phase relation usually deviates from in-phase 
properties when the disturbances are reflected off of walls or exhibit 
standing -mode properties of acoustic resonance. Such disturbances should 
grow less rapidly than traveling wave disturbances. 

CONCLUDING REMARKS 

This analysis has shown the response properties of a convectively con- 
trolled burning process to be highly sensitive to the harmonic distortion 
of a flow field disturbance. Harmonic distortion can amplify or attenuate 
the response in a manner which should significantly alter the growth or 
decay of a disturbance. Directing more attention toward these response 
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properties of convectively controlled turning could increase the understand- 
ing of the dynamic behavior of many reactive systems and provide new ap- 
proaches to experimentally modify the behavior of a specific system. 
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APPENDIX - SYMBOLS 

speed of sound 
characteristic dimension 
harmonic order 



magnitude of relative velocity vector 
gas velocity 

harmonic coefficient for velocity disturbance of order n 
steady relative velocity of particle or surface with respect to 
the gas 

dimensionless energy release rate 

energy release rate per unit time and volume 

ratio of specific heats • 

harmonic phase angle for velocity- pres sure relation 
gas density 

harmonic coefficient for density disturbance of order n 
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cp n harmonic phase angle 

4 gas viscosity 

to frequency 

' primed quantities denote dimensionless perturbations, x' = x - x/ 


barred quantities denote mean values 
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Figure 4. -Effect of second harmonic content on response and wave 
shape. q >2 s 0, = 83 = 0, Av/c * 0. 02, y = 1.2. 
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Figure 5. - Response and wave shapes for muttlharmonlc 
disturbances. <p ■ 0, 6 - 0, y * 1. 2. 
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Figure 6. - Effect of harmonic phase angle on the response to 
multfharmonic disturbances equation (9). 0*0, Av/c = 
0.02, y* 1.2. 
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